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Half-life estimates for neutrinoless double beta decay depend on particle physics models for lepton- 
flavor violation, as well as on nuclear physics models for the structure and transitions of candidate 
nuclei. Different models considered in the literature can be contrasted — via prospective data — with 
a "standard" scenario characterized by light Majorana neutrino exchange and by the quasiparticle 
random phase approximation, for which the theoretical covariance matrix has been recently esti- 
mated. We show that, assuming future half-life data in four promising nuclei (^''Ge, ^^Se, "°Te, and 
^\ • ^'^''Xe), the standard scenario can be distinguished from a few nonstandard physics models, while 

' being compatible with alternative state-of-the-art nuclear calculations (at 95% C.L.). Future signals 

, in different nuclei may thus help to discriminate at least some decay mechanisms, without being 

' spoiled by current nuclear uncertainties. Prospects for possible improvements are also discussed. 

^ PACS numbers: 23.40.-s, 02.70.Rr, ll.30.Fs, 12.60.-i 

\0 [ I. INTRODUCTION 

. The search for lepton number violation has gained new momentum after the discovery of neutrino flavor transfor- 
CIh' mations driven by ly masses and mixing [Tj. In particular, the search for neutrinoless double beta (Oi^/?/3) decay. 
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which violates the lepton number by two units, is at the forefront of current neutrino physics 0. 
For a given candidate nucleus i = (Z, A), the 01/(3(3 half- life Ti can be expressed as 

>: t-'^gUm^{'>4, (2) 



where is an accurately calculable kinematical phase-space factor, M[ is the 0i^/3/9 nuclear matrix element (NME), 
embedding the nuclear physics aspects of the process, and Xi is a lepton number violation parameter, embedding the 
particle physics aspects of the process. The index £ labels different possible mechanisms underlying 0iyf3(3 decay 
' If the three known neutrinos i^i are described by Majorana (rather than Dirac) 4-spinors, then their masses and 
the Ve mixing matrix elements ILi can generate one such mechanism, with Xi being equal to the so-called "effective 
Majorana neutrino mass," m^^ 



1=1 



(3) 



Apart from the above "standard" mechanism with the exchange of light Majorana neutrinos, Qvj3(3 decays can also 
be generated by new particles and interactions (and even new dimensions) beyond those of the standard electroweak 
model, as recently reviewed in [1, [^. 

Different models of particle and nuclear physics affect half-life predictions Ti in different ways, and may thus be 
discriminated, in principle, if the overall uncertainties are smaller than the spread of the predictions. In particular, 
the use of half- life ratios, 

^ G] \M^:\' 

Tj Gi\Mf\^ ' ^ ' 

has been advocated in order to effectively reduce the dependence on the following: (i) unknown particle physics 
parameters via cancellation of the factor A^; and (ii) nuclear physics uncertainties 0, via cancellation of 

systematics in the NME ratio, so as to enhance the discrimination power of future data. In this context, the work Q 
represents (to our knowledge) the only quantitative study of particle and nuclear physics model discrimination, based 
on prospective multi-isotope data and on guessed uncertainties — usually dominated by theoretical NME errors. 

Recently, there has been significant progress toward a reduction and a better understanding of the theoretical 
nuclear uncertainties. In particular, state-of-the-art NME calculations based on the quasiparticle random phase 
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TABLE I: Prospective half-life sensitivities at 90% C.L. {Tf°) for different nuclei i in promising future projects, as reported in 
All projects plan a second phase with lower backgrounds and higher sensitivies. 



i 




Project 




2.0 X 10^*^ 


GERDA, MAJORANA 


«^Se 


2.0 X 10^® 


SuperNEMO 




2.f X fO^^ 


CUORE 


"6xe 


6.4 X 10^^ 


EXO 



approximation (QRPA) [71, |8| or on the Shell Model [9| are increasingly converging within their estimated errors, at 
least in the most studied Oi/(3P case with light Majorana neutrino exchange. Moreover, within a specific QRPA model 
and its variants [2], the theoretical error matrix has been estimated in detail 10], showing that high — and previously 
ignored — correlations pij exist between NME errors in any pair of nuclei (i, j). High correlations imply that NME 
errors are largely cancelled in ratios Ti/Tj, thus motivating a novel statistical analysis of prospective multi-isotope 
data where the pij are explicitly included — which is the goal of this work. 

The structure of our paper is as follows. In Sec. II we describe in detail our "standard scenario," characterized 
by detailed QRPA estimates of NME's for light Majorana neutrino exchange [13] ■ This scenario is used as null 
hypothesis in statistical tests, based on prospective data for four Oi^/3/3 candidate nuclei (^^Ge, ®^Se, ^^°Te, and 
^^^Xe). In Sec. Ill we assume mock data according to alternative particle physics model expectations, and show that 
two (out of seven) models can reject the null hypothesis at > 95% C.L. In Sec. IV we repeat the exercise with different 
nuclear physics models, and show that none of them can reject the null hypothesis at 95% C.L. Therefore, future 
0^/3/3 signals may indeed start to discriminate some underlying particle physics mechanisms, without being spoiled 
by currently estimated nuclear physics uncertainties. Theoretical error correlations are proven to be crucial in this 
context. Conclusions and prospects for future improvements are discussed in Sec. V. 

II. NULL HYPOTHESIS AND STATISTICAL APPROACH 

In this work, we focus on four candidate nuclei for Qv[3l3 decay, which will be probed in promising, next-generation 
experiments IT]: ''^Ge, *^Se, ^^°Te, and ^'^^Xe. Table I shows the expected half-life sensitivity in the first experimental 
phase, which might be followed by an even more sensitive second phase [lH. Among these four nuclei, ''^Ge is perhaps 
the most studied, both theoretically and experimentally. Therefore, its half-life Tj may be used as an appropriate 
benchmark 0, namely, as a common denominator in ratios Ti/Tj. For definiteness, we assume a prospective value 

T(^^Ge) = 10^6 y , (5) 

reminding that such absolute half-life is immaterial and provides just a normalization. Our main results do not depend 
on this specific choice for the "denominator" nucleus and its half-life. It is also appropriate to linearize the analysis 
by taking logarithms of half- lives [13] , 

T, =logio(T,/y) , (6) 

so as to deal with differences (t^ — Tj) rather than with ratios {Ti/Tj). 

Recently, state-of-the-art calculations of Qv(3j3 NME (within the QRPA and for light Majorana neutrino exchange 
) have been thoroughly investigated from a statistical viewpoint [1^, in order to assess their inherent theoretical 
uncertainties. In particular, it has been shown that the NME calculations display not only large variances (as is 
well known), but also large covariances among any two nuclei, which can be accounted for by a correlation matrix 
Pij. The errors and correlations estimated in [lO| originate from the following variants in the QRPA ingredients 
used as inputs: {i) two values for the axial coupling, namely, gA = 1-25 (bare) and qa = 1-00 (quenched); {ii) two 
approaches to short-range correlations (s.r.c), the so-called Jastrow-type s.r.c, and the unitary correlation operator 
method (UCOM); (iii) three sizes for the model basis, small, intermediate, and large; {iv) two many-body models, 
namely, QRPA and its renormalized version. All the 24 variants include the NME uncertainties induced by fitting 
the particle-particle strength parameter g^p via the experimentally observed 2vj3(3 decay process [3] • 

The final results of [To^, together with the normalization in Eq. (O, can be translated easily into a set of predictions 
for the logarithmic half-lives and their ±1ct errors. 



r° ± si (with correlations pij) , 



(7) 



3 



TABLE II: Null hypothesis: central values = logj^Q(ri/y), together with their theoretical errors and correlations pij, in our 
standard QRPA scenario with light Majorana u exchange 10], normalized to a benchmark half-life Ti/y — 10^^ for i = ^'^Ge. 



Null hypothesis 


Ge 


82 o 

Se 


Tc 


136 V 

Xe 


QRPA, standard rf 


26.000 


25.480 


25.430 


25.888 


1-cr error 


0.244 


0.270 


0.316 


0.374 


Correlations pij 


1.000 


0.978 


0.899 


0.805 






1.000 


0.927 


0.846 








1.000 


0.916 



1.000 



where the superscript (°) is intended to mark the "null hypothesis" in the following statistical tests. Table II shows 
the resulting numerical values for our null hypothesis (or "QRPA standard scenario"). 

The null hypothesis can be tested via prospective experimental data. We assume that several experiments will 
observe positive Qv(3(3 signals, and determine the decay lifetimes with uncorrelated experimental errors, 

Ti ± Si (data) . (8) 

Next-generation experiments might aim at a half- life accuracy of STi/Ti ~ 20% j2|], namely, 

s,; ~ 0.08 , (9) 

which we adopt hereafter. The specific value of Si is not crucial, as far as theoretical uncertainties remain dominant 
(s^ ;3> Si). If the differences between standard predictions and mock data (rf — Ti) are significantly larger than the 
theoretical errors (s^) in at least one nucleus (other than ^^Ge, where — Ti — 26.000 by construction), the null 
hypothesis will be rejected — or it will be accepted otherwise. 
The differences can be properly evaluated by statistics [loj . 

N 

X'^ 5](r.-r°)I^.,(r,-r°), (10) 
where the inverse error matrix W includes both theoretical errors (correlated) and experimental errors (uncorrelated) : 

m^'^P^,s",s",+5.,S,S, . (11) 

In the above expression, iV is the number of nuclei considered in the analysis. We shall take = 4 by default, but 
also will show variations for N = 3, whereas the case with N = 2 nuclei (providing a single half-life ratio with relatively 
little constraining power) will not be commented in this work. Since the ^^Ge half-life is used for normalization in all 
cases, the correct number of degrees of freedom is — 1. The rejection probability associated to the test [l^ 
is thus 

Pr=Pr{x^;N~l) . (12) 

Moreover, we set a threshold Pr > 95% for a statistically significant rejection of the null hypothesis (implying, e.g., 
> 7.8 for A^ = 4). 

In comparison with the statistical approach of ^5] , our method uses a simple x^ statistics (rather than Monte Carlo 
simulations), and includes theoretical error correlations (which will be proven to be crucial). On the other hand, we 
test a single null hypothesis (the standard QRPA case, where such correlations are well defined [13 )i and focus on a 
smaller set of candidate nuclei than in 5]. In this sense, both our results and those in provide quantitative, but 
complementary, assessments of future tests of models for Oi/P/3 decay. Prospects for more comprehensive analyses, 
including the merits of both approaches, are discussed below in Sec. V. 



III. ANALYSIS OF PARTICLE PHYSICS MODELS 



Neutrinoless double beta decay can be triggered by several lepton- flavor mechanisms, involving either long-range 
or short-range interactions. Although two or more mechanisms might be present (and interfere) at the same time, 
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TABLE III: Estimated values of Ti = logiQ(Ti/y) in various particle physics models for Qvj3p decay (different from light 
Majorana v exchange), normalized to a benchmark half-life Ti/y — for i = ^''Ge. 



Particle physics model 


™Gc 




»"Se 




"«Xe 


Refs. 


Heavy v 


26.000 




25.415 


25.146 


25.591 


fl41 


SUSYtt 


26.000 




25.431 


25.462 


25.863 




SUSY g 


26.000 




25.447 


25.230 


25.724 


[16] 


RHC n 


26.000 




25.462 


25.301 


25.732 


[17] 


RHC A 


26.000 




25.146 


25.114 


25.826 


[17] 


KK+1 


26.000 




25.380 


25.279 


26.519 


[18] 


KK-1 


26.000 




25.415 


25.255 


25.892 


fl8] 


TABLE IV: Test of the null hypothesis (QRPA, standard), assuming Qv(3l3 half- life data as predicted by different particle 
physics models, with 20% experimental uncertainties. Second column: rejection probability Pr (%) with four nuclei. Third 
column: rejection probability range for three nuclei (^^Ge plus any other two). The null hypothesis can be rejected at > 95% 
C.L. in two models, marked by Fourth and fifth columns: as in the previous two columns, but without correlations. 


Particle physics model 


4 nuclei, % 


3 nuclei, % 




4 nuclei, % 


3 nuclei, % 












(no pij) 


(no pij) 


Heavy v 




53.3 


49.6 - 71.3 




29.8 


28.0 - 49.4 


SUSY TT 




5.2 


5.6 - 14.7 




0.2 


0.7 - 2.0 


SUSY ~g 




27.1 


18.8 - 47.7 




9.8 


9.4 - 24.4 


RHC 77 




10.5 


17.9 - 24.1 




4.5 


7.7 - 14.9 


RHC A * 




97.6 


89.5 - 97.9 




50.1 


38.3 - 69.1 


KK+1 * 




99.9 


35.4 - 99.9 




61.8 


15.7 - 77.0 


KK-1 




37.0 


14.9 - 57.2 




4.8 


2.6 - 15.7 



we shall restrict ourselves to cases with just one dominant contribution. In this context, several models have been 
recently reviewed in Refs. 0,01, to which we refer the reader for details and references to earlier literature We 
focus on the following nonstandard mechanisms for OvPf] decay: 

(1) Exchange of Majorana neutrinos with heavy masses (m^ > 1 GeV) (Heavy z/); 

(2) Supcrsymmctric models with R-parity violation and pion exchange [l5j (SUSY tt); 

(3) As above, but with gluino exchange [T^ (SUSY g); 

(4) Left-right symmetric models with leptonic right-handed currents (RHC) coupled to hadronic left-handed currents 
[13 (RHC r?); 

(5) As above, but coupled to hadronic right-handed currents [T7| (RHC A); 

(6) Kaluza-Klein neutrino exchange via extra dimension with radius R = (300 eV)^^ and brane-shift parameter (in 
GeV-i): a = 10+1 (KK+l); 

(7) As above, but with a = IQ-^ (KK-1). 

Another SUSY-accompanied mechanism listed in [3| provides basically the same predictions as the SUSY g model in 
our four reference nuclei, and thus it is not separately considered here. 

Table HI shows the estimated logarithmic half-lives in the above particle physics models, as derived from the 
numerical compilations in 0, [H], but with the normalization in Eq. ([5]). For each model, we assume the r^'s as central 
values of mock data, with a putative experimental uncertainty Si as in Eq. We then test the null hypothesis in 
Table II via the test in Eq. 
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Before discussing quantitative results, we observe that the spread induced by different models in Table III is 
often comparable or smaller than the theoretical uncertainties s^. However, there are also some outliers, e.g., a very 
high ^^^Xe half-life prediction in the KK+1 model. Therefore, it should be possible to distinguish, at least, a few 
nonstandard physics models from the standard case of light Majorana neutrino exchange (represented by the null 
hypothesis). 

Table IV shows the results of our statistical analysis. Using all four nuclei, the rejection probability for the 
null hypothesis ranges from an insignificant 5.2% (SUSY tt model) to a highly significant 99.9% (KK+1 model). 
Also another model (RHC A) provides a relatively high Pr — 97.6%. Therefore, with four nuclei, at least these two 
nonstandard 0i>P(3 mechanisms can be distinguished from the standard one at > 95% C.L., while the other mechanisms 
are statistically indistinguishable from the null hypothesis. 

By discarding one nucleus (except the benchmark ''^Ge one), the rejection level for — 1 = 2 degrees of freedom 
may change considerably in either directions, depending on the role of that nucleus in the total x^- In general, it 
may be advantageous to drop a nucleus whose standard and nonstandard half-life predictions are similar. Viceversa, 
dropping a nucleus with rather different predictions may weaken the statistical power of the test. 

Table IV also shows the range of rejection probabilities Pr with three nuclei. It turns out that the RHC A and 
KK-l-1 models can reach P > 95% not only with four nuclei, but even with some three-nuclei combinations. In 
particular, we find that 95% C.L. can be exceeded by combining the benchmark Ge datum with: (i) either Se-|-Te 
or Se-|-Xe data for the RHC-A model, and {ii) either Xe-|-Se or Xe-|-Te data for the KK+1 model. In any case, with 
either three or four nuclei, correlations of theoretical errors are crucial for the statistical power of the test: if they 
were neglected {pij = Sij ) , no model could be distinguished from the null hypothesis at a statistically significant level 
(95% or higher), as reported in Table IV. 

The role of theoretical error correlations can be better appreciated by means of Fig. 1, which charts the logarithmic 
half-lives for all pairs of nuclei. The slanted ellipses represent the null hypothesis (QRPA, standard), while the 
crosses represent the assumed mock data for the two most deviant models (RHC A and KK+1), within la errors. 
By construction, all central values for ^^Ge are aligned at T = 10^^ y (vertical dotted lines). The relatively large 
"distance" between "theory" (ellipses) and "data" (crosses) depends crucially on the ellipse errors and orientation. 
If correlations were dropped (pij = Sij), all the ellipses would cover a much larger fraction of the planes, and their 
frontier would get closer to the crosses, in general. Viceversa, if correlations were maximal {pij = 1), the ellipses 
would shrink to slanted segments: theoretical uncertainties then would exactly cancel in half-life ratios, and the 
analysis would be dominated by experimental errors. This fact underlines the importance of detailed estimates of the 
theoretical covariance matrix, in addition to other motivations discussed in p^ . 

In conclusion, it appears that future measurements of 0v(3P decay half-lives in four (or even three) nuclei can 
discriminate the standard mechanism of light Majorana neutrino exchange from at least the most deviant nonstandard 
mechanisms. The theoretical covariance matrix plays a crucial role in the statistical analysis, and its control might 
improve the discrimination power of future measurements. 



The particle physics models considered in the previous Section (Tab. Ill and IV) have been implemented, in 
the available literature, within different approaches to the nuclear physics structure and transitions for the QvPfi 
candidate nuclei considered. Therefore, it makes sense also to analyze the spread of predictions due to various nuclear 
approximations. In this Section, we focus on the standard Qu(3(3 physics mechanism (light Majorana v exchange), and 
on a set of recent, state-of-the-art calculations of the associated nuclear matrix elements: 

(1) QRPA calculations by Suhonen and Kortelainen (S&K) for fixed qa = 1-00 and Jastrow s.r.c; 



(6) SheU Model with gA 1.25 and UCOM s.r.c. 

See also the discussion in the Appendix of [l3], where the NME from the above nuclear models are compiled and 
compared with the NME used in our null hypothesis. We mention here that a third approach, based on the microscopic 
interacting boson model (IBM) T^, provides NME's in agreement with QRPA; however, the recent results in do 
not include ^^^Xe estimates and thus are not included in this work. 



IV. ANALYSIS OF NUCLEAR PHYSICS MODELS 



(2) 
(3) 

(4) 
(5) 



QRPA (S&K) for qa = 1.00 and UCOM s.r.c; 
QRPA (S&K) for gA = 1.25 and Jastrow s.r.c; 
QRPA (S&K) for gA = 1.25 and UCOM s.r.c; 



Shell Model calculations [9j for fixed gA — 1-25 and Jastrow s.r.c; 
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T(-Ge)/y T(-Se)/y T(-°Te)/y 



FIG. 1: Projections of logarithmic half-life estimates in the coordinates planes of each pair of nuclei. The ellipses represent 
our null hypothesis (QRPA, standard), with evident correlations of theoretical errors (at la). The crosses represent mock data 
(with assumed 20% accuracy at la) centered on the predictions of two nonstandard particle physics models described in the 
text, namely, RHC A and KK+l. Data and predictions are conventionally aligned at a benchmark value T = 10^^ y for ^^Ge 
(vertical dotted lines) 
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TABLE V; Estimated values of n — logjQ(ri/y) in state-of-the-art nuclear physics models for GvPfS decay (for light Majorana 
V exchange), normalized to a benchmark half-life Ti/y — for i = ^^Ge. 



Nuclear physics model 


^«Ge 


«2Se 


130 Te 


i^'^Xe 


Refs. 


QRPA (S&K), gA = 1.00, Jastrow 


26.000 


25.669 


25.351 


25.625 


[8] 


QRPA (S&K), QA = 1.00, UCOM 


26.000 


25.663 


25.319 


25.621 


m 


QRPA (S&K), QA = 1.25, Jastrow 


26.000 


25.679 


25.415 


25.713 


m 


QRPA (S&K), QA = 1.25, UCOM 


26.000 


25.669 


25.363 


25.689 




Shell Model, qa = 1.25, Jastrow 


26.000 


25.399 


25.227 


25.359 


[9] 


Shell Model, qa = 1.25, UCOM 


26.000 


25.407 


25.207 


25.343 




TABLE VL Test of the null hypothesis (QRPA, standard), assuming QufSp half-life data as predicted by different nuclear physics 
models, with 20% experimental uncertainties. Second column: rejection probability Pr (%) with four nuclei. Third column: 
rejection probability range for three nuclei (^^Ge plus any other two). The null hypothesis is not rejected, within 95% C.L. 
Fourth and fifth columns: as in the previous two columns, but without correlations of theoretical errors. 


Nuclear physics model 


4 nuclei, % 


3 nuclei, % 




4 nuclei, % 

{P^J = 0) 


3 nuclei, % 
(p,:, = 0) 


QRPA (S&K), QA = 1.00, Jastrow 


86.8 


45.0 - 93.7 




19.4 


22.5 - 37.0 


QRPA (S&K), QA = 1.00, UCOM 


87.0 


43.0 - 93.1 




20.5 


23.6 - 36.5 


QRPA (S&K), gA = 1.25, Jastrow 


79.9 


30.8 - 90.1 




12.9 


10.0 - 29.8 


QRPA (S&K), QA = 1.25, UCOM 


80.9 


28.1 - 89.9 




14.2 


14.5 - 30.3 


Shell Model, qa = 1.25, Jastrow 


78.9 


46.4 - 89.4 




50.3 


21.0 - 68.4 


Shell Model, qa = 1.25, UCOM 


81.1 


52.5 - 90.7 




53.6 


23.5 - 71.3 



Table V shows the estimated logarithmic half-lives Ti in the above six nuclear models, as derived from the compilation 
in [To| via the normalization in Eq. ([5]). Analogously to the previous Section, for each model we assume the r^'s as 
central values of mock data with errors st [Eq. ©], and perform a test of the null hypothesis of Table II. 

Table VI shows the results of the statistical test, with either four or three nuclei, with and without theoretical 
error correlations. It appears that, in any case, none of the six nuclear models can reject the null hypothesis (i.e., 
our reference QRPA calculations [3]) at > 95% C.L. In other words, at such level of significance, the six nuclear 
models considered are phenomenologically indistinguishable from our reference QRPA model. This result is not in 
contradiction with the significant differences (up to ^ ia) existing among the corresponding NME's [10|]. In fact, 
we are not really dealing with the absolute, unobservable NME, but only with their relative sizes (with respect to 
the benchmark ^^Ge expectations), which exhibit a smaller spread. Of course, the larger spread of absolute NME 
will reappear elsewhere, namely, as a nuclear model uncertainty in the lepton flavor violation parameter A^, whose 
reconstruction uncertainties are beyond the scope of this work (see for a discussion). 

In Table VI, the two independent models which deviate most from the null hypothesis are: QRPA (S&K) model 
with gA = I and Shell Model with gA = 1.25, both with UCOM s.r.c. Figure 2 shows mock data according to such 
models (crosses), as compared with our standard QRPA expectations (ellipses) in the same coordinate planes as in 
Fig. 1. The comparison of Figs. 1 and 2 shows at a glance that: (?) the two most deviant nuclear physics model 
(Fig. 2) are closer to standard predictions than the two most deviant particle physics models (Fig. 1); and (ii) in 
any give panel, the largest deviations between theory and data are generally opposite in Figs. 1 and 2, implying that 
nuclear physics variations do not mimic particle physics variations in this context. 

We can thus conclude that, using prospective OvP/S data in four (and possibly just three) promising nuclei, the 
discrimination of some particle physics mechanisms (in particular, RHC A and KK-t-1 in our analysis) is not spoiled 
by estimated uncertainties in state-of-the-art nuclear theory calculations. This is the main result of our work. 



^ QRPA, standard 
QRPA (S & K) 
Shell Model 




^q27 ^q25 



T(-Ge)/y 



T(-Se)/y 



10'° 10 
T("°Te)/y 



27 



FIG. 2: As in Fig. 1, but with mock data centered on the predictions (for light Majorana v exchange) of two alternative nuclear 
models, namely: QRPA (S&K) with gA ~ 1, and Shell Model with gA ~ 1.25, both with UCOM s.r.c. See the text for details. 
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V. SUMMARY AND DISCUSSION 

In this paper, we have focused on a reference Qv[3l3 scenario characterized by hght Majorana neutrino exchange 
(concerning particle physics) and by QRPA (concerning nucfear physics), for which the theoretical covariance matrix 
can be defined in detail We have used this scenario as a null hypothesis to be tested with prospective half-life 
data (with 20% accuracy) in four promising nuclei (™Ge, ^^Se, -^^°Te, and -^^^Xe), with the first nucleus being used as 
a conventional benchmark. Through our statistical approach, one can properly quantify and implement the correlated 
uncertainties associated to half-life ratio tests of decay, which were proposed, e.g., in [3,0] for particle physics 

models and in [^Q for nuclear physics models. 

By setting a 95% C.L. threshold for statistical significance, we have shown that two nonstandard models for 
neutrinoless double beta decay, based on right-handed leptonic currents and on Kaluza-Klein excitations in extra 
dimensions, can reject the null hypothesis. Conversely, state-of-the-art nuclear physics models are consistent with the 
null hypothesis. Therefore, in case of positive 0vP(3 signals, next-generation experiments will have some discrimination 
power about the underlying particle physics mechanism, despite known nuclear physics uncertainties. Since theoretical 
uncertainties dominate the analysis, their reduction can only increase the discrimination power of future Qv(3[3 data. 
Correlations among theoretical errors appear to be crucial ingredients in this context. 

Our approach can be improved and extended in several ways, which we list in order of increasing difhculty. Further 
prospective data can be easily included: the analysis does not need to be limited to four nuclei (see, e.g., 5,]) and, 
in case of negative results, it might include one-sided limits besides two-sided signals. Analogously, further particle 
or nuclear physics models can be included and tested. However, improvements on model discrimination power are 
less obvious, since theoretical covariances have been evaluated only for a single set of QRPA calculations and for 
light Majorana neutrino exchange ^IQ] . For instance, in order to compare any two models as proposed in [5] (rather 
than each model to the same null hypothesis), one should properly attach a covariance matrix to each model. To 
reach this goal, future NME estimates should include a "statistical set" of variants due to different input choices in 
their model parameter space, as also emphasized in [l^]- It would also be useful to systematically evaluate NME's 
for various particle physics mechanisms in one and the same nuclear physics model [20| (and vice versa) in order 
to gauge more precisely the impact of each model on observable half-lives, and to investigate how many nuclei are 
needed for model discrimination at a given C.L. Needless to say, these goals imply a long-term theoretical research 
program, which will be probably fully pursued only after unmistakable signals for Qv(3(3 decay are found in two or 
more nuclei. It is encouraging, however, to observe that, at least for the standard case of light Majorana neutrino 
exchange, very different approaches (QRPA 0, M: Hi, Shell Model Q and IBM [H) are converging more than it 
could be hoped just a few years a go. Furthermore, within the QRPA, the convergence improves as further nuclear 
structure constraints are included j2ll. [2^. First-order electroweak reactions, probing the so-called first and second 
leg of QvPP decay, might also help to reduce parametric uncertainties of nuclear models [13, HI] . 

Finally, when {)v(3(3 half-life data will be hopefully established in multiple nuclei, additional tests of the underlying 
particle physics mechanism might include further constraints from the kinematical distributions of the emitted elec- 
trons [3, [25[ , from the comparison of decay with electron capture jl^ , from branching ratios of Qvfif} decays 
to g round and excited states [l^], and from possible hnks with other lepton- flavor violating processes (e.g., ^ — > e) 
[28| . Such a rich phenomenology will then warrant more comprehensive analyses than the one proposed in this work. 
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